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Cl Abstract 

This article presents the convergence analysis of a sequence of piecewise 
constant and piecewise linear functions obtained by the Rothe method 

1 i to the solution of the first order evolution partial differential inclusion 

u'(t) + Au(t) + l* d.J{iu(t)) B f(t), where the multivalued term is given by 
O^l the Clarke subdifferential of a locally Lipschitz functional. The method 

provides the proof of existence of solutions alternative to the ones known 
On in literature and together with any method for underlying elliptic prob- 

lem, can serve as the effective tool to approximate the solution numeri- 
cally. Presented approach puts into the unified framework known results 
for multivalued nonmonotone source term and boundary conditions, and 
generalizes them to the case where the multivalued term is defined on 
the arbitrary reflexive Banach space as long as appropriate conditions are 
satisfied. In addition the results on improved convergence as well as the 
numerical examples are presented. 

> 

k*n 1 Introduction 
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Partial differential inclusions with the multivalued term given in the form of 
Clarke subdifferential are known as hemivariational inequalities (HVIs). HVIs 
are the natural generalization of the inclusions with monotone multivalued term 
(which lead to variational inequalities) and were firstly considered by Pana- 
giotopoulos in early 1980s. For the description of the origins of HVIs and un- 
derlying mathematical theory we refer the reader to the book [25] . 
This paper deals with the first order evolution inclusion of type u'(t) + A(u(t)) + 
i*dJ{bu(t)) 9 f(t). Such problems are known as parabolic HVIs or boundary 
parabolic HVIs depending whether an operator t is the embedding operator 
from H 1 ^) to L 2 (tt) or the trace operator from H 1 ^) to H^(dn). The first 
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case corresponds to multivalued and nonmonotone source term in the equation 
and the second one to multivalued and nonmonotone boundary conditions of 
Neumann-Robin type. Such inclusions are used to model the diffusive transport 
through semipermeable membranes where the multivalued term represents the 
semipermeability relation [24] and the temperature control problems where the 
multivalued term represents the feedback control [IB], [15j . 
The existence of solutions to problems governed by inclusions of considered type 
was investigated by many authors. There are several techniques used to obtain 
the existence results: 

• Classical Faedo-Galerkin approach combined with the regularization of 
the multivalued term by means of a standard mollifier; solutions of un- 
derlying system of ordinary differential equations are proved to converge 
(in appropriate sense) to the function which is shown to be the solution 
of analyzed HVI. This technique was used in context of parabolic HVIs 
by Miettinen [22] , Miettinen and Panagiotopoulos [23] and Goeleven et al. 

[HJ. 

• The approach based on the notion of upper and lower solutions. The 
solution is shown to be the limit of solutions of problems governed by the 
equations obtained by the regularization of the multivalued term together 
with the truncation by the lower and upper solutions. The distinctive 
feature of this approach is that the growth conditions on the multivalued 
term are replaced by the assumption of the existence of lower and upper 
solutions. The technique was used for parabolic HVIs by Carl [3] and 
developed in 0, [BJ, [7], 0. 

• The technique based on showing that the analyzed HVI satisfies the as- 
sumptions of the general framework for which the appropriate surjectivity 
result holds. This approach was used by Liu [19] and by Migorski |26| and 
developed for the boundary case in |27j . 

• The technique based on adding to the inclusion the regularizing term 
multiplied by e > 0, showing that the solutions to obtained problems 
satisfy some bounds uniformly in e and passing to the limit e — > 0. This 
technique was used for parabolic HVIs by Liu and Zhang [IT] and Liu [IB] 
and developed in [20 ] . [2T ] , 

It should be remarked that above techniques are either nonconstructive (i.e. 
they are based on surjectivity result) or constructive but not effective (i.e. re- 
quire a priori knowledge of lower and upper solutions, or require additional or 
smoothing terms in the problem). 

In contrast to the existence theory, numerical methods to approximate effec- 
tively the solutions to parabolic HVIs were not considered by many authors. 
In the book of Haslinger, Miettinen and Panagiotopoulos [IBJ the convergence 
of solutions obtained by the finite element approximation of the space variable 
and finite difference approximation of the time variable is proved. However 
only the case of the linear operator A and the multivalued source term (and 
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not boundary conditions) is considered (see Remark 4.10 in [IS]). In [15] the 
authors proved the convergence of the finite difference scheme (with respect 
to both time and space variable) for the case of multivalued source term (i.e. 
U = H in the sequel). 

Our approach uses the so-called Rothe method (known also as time approxi- 
mation method) and allows to extend any numerical method that is used to 
solve the stationary, elliptic inclusions with the multivalued term given as the 
Clarke subdifferential, to time dependent, parabolic problems. The key idea 
is the replacement of time derivative with the backward difference scheme and 
solve the associated elliptic problem in every time step to find the solution in 
the consecutive points of the time mesh. It is proved that the results obtained 
by such approach approximate the solution of the original problem. 
On the other hand, the Rothe metod provides the proof of existence of solutions. 
In contrast to other approaches this metod, as long as one can solve underlying 
elliptic problems, does not require any smoothing or other additional regulariz- 
ing terms in the inclusion. Furthermore the presented approach allows to study 
the inclusions with multivalued term given on the domain and on the domain 
boundary within the unified framework in which the multifunction that appears 
in the problem is defined on an arbitrary reflexive Banach space, which satisfies 
the appropriate assumption {H{U) in the sequel). This assumption is proved 
to generalize the case of inclusions with multivalued boundary conditions and 
the ones with multivalued source term (see Section [3] and examples of problem 
settings in Section [8]). 

The Rothe method for parabolic nonlinear PDEs with pseudomonotone opera- 
tors is described in the monograph of Roubicek [33] , where also the results for 
the monotone multivalued problems are presented. In the context of parabolic 
HVIs the variant of the Rothe method was used to show existence of solutions 
to problems with hysteresis in [33] and [25] , but there only the case of linear op- 
erator A and / € L 2 (0, T; L 2 {Vt)) (which excludes nonhomogeneous Neumann 
conditions) was considered and besides only the case of the multivalued and 
nonmonotone term source term was analyzed. 

In Section [2] some basic definitions are recalled. Section [3] presents the gen- 
eralization of the Lions- Aubin Compactness Lemma that justifies the usage of 
the assumption H(U) in the sequel. Problem setup and the assumptions are 
presented in Section [4] The auxiliary elliptic problems solved in every time 
step, which are the key idea of the Rothe method, are formulated and analyzed 
in Section [5] Convergence of piecewise linear and piecewise constant functions 
constructed basing on the solutions of auxiliary problems as well as the fact that 
the limit solves the original problem is proved in Section [6 Some stronger con- 
vergence and uniqueness results are established in Section 7\ Finally in Section 
[8] it is shown that the cases of multivalued boundary condition and source term 
are the special cases of presented general framework and a simple numerical 
example is delivered. 
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2 Preliminaries 

In this section we recall several key definitions that will be used in the sequel. 
For a locally Lipschitz functional j : X K, where X is a Banach space, gen- 
eralized directional derivative (in the sense of Clarke) at x £ X in the direction 
z E X is defined as 

f(x;z) = hmsup feLMzM . 

2/->-x,A->0+ A 

Generalized gradient of j (in the sense of Clarke) is the multifunction dj : X — > 
2 X * defined by 

0j(aO = G X* : j°(x;y) > for all y e X}, 

where (■, ■) stands for the duality pairing between X and X* . For the properties 
and the calculus of the Clarke gradient see [TU] . 

Recall that the multifunction A : X —> 2 X , where X is a real and reflexive 
Banach space is pseudomonotone if 

(i) A has values which are nonempty, weakly compact and convex, 

(ii) A is use from every finite dimensional subsepace of X into X* furnished 
with weak topology, 

(iii) if v n — > v weakly in X and u* e A(v n ) is such that limsup J1 _ i . 00 (w* , v n — 
v) < then for every y € X there exists u(y) G ^4(w) such that (u(y), v — 
y) < liminfn^oo (v*,v n - y). 

Note that sometimes it is useful to check the pseudomonotonicity of an op- 
erator via the following sufficient condition (see Proposition 1.3.66 in or 
Proposition 3.1 in [8]). 

Proposition 1. Let X be a real reflexive Banach space, and assume that A : 
X — > 2 X satisfies the following conditions 

(i) for each v € X we have that A(v) is a nonempty, closed and convex subset 
ofX*. 

(ii) A is bounded. 

(iii) If v n — ¥ v weakly in X and w* — > v* weakly in X* with v* € A(v n ) and if 
lim sup Tl _j. 00 (u*, v n — v) < 0, then v* € A(v) and (w* , v n ) — > (v* , v) . 

Then the operator A is pseudomonotone. 

We also recall (see for instance Proposition 1.3.68 jT^]) that the sum of two 
pseudomonotone multifunctions is pseudomonotone. 
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3 Generalization of Lions- Aubin Lemma 

For a Banach space X, 1 < p < oo and a finite time interval I = (0, T) we 
consider the standard spaces L P (I;X). Furthermore we denote by BV(I;X) 
the space of functions of bounded total variation on /. Let it denote any finite 
partition of 7 by a family of disjoint subintervals {<Tj = (o-i,bi)} such that 
/ = 1J" =1 <Ji- Let T denote the family of all such partitions. Then we define the 
total variation as 



IMI w(/ : x) = sup \ V ||x(6i) - x{ai)\\ x ) ■ 
As a generalization of above definition for 1 < q < 00 we can define a seminorm 



\x\\ q 

< "BVi(I;X) 



= sup <^ \\ x ( b i) - x ( a i)\\x \ ■ 



For Banach spaces X 1 Z such that IcZwe introduce a vector space 

M p ^(I; X, Z) = L p (I; X) n BV q (I; Z). 

Then M p ' q (7; X, Z) is also a Banach space for 1 < p, q < 00 with the norm 
given by || • \\ L p(i;X) + II ■ \\bv(I-,z)- 

Let us recall Theorem 1 of [33] (see also Theorem 1 of 32J and Proposition 2.1 

of El]). 

Theorem 1. Let 1 < p < 00 and X be a real Banach space. A subset Q C 

L p (0,T; X) is relatively compact in a Banach space L p (0,T; X) provided the 
following two conditions hold 

• for every < t% < t<i < T the set 

ft 



G{t l: t 2 ) ■■= jjf u(t)dt:ueG^ 



is relatively compact in X , 
• Q is strongly integrally equicontinuous i.e. 

limsup/ \\u{t + h) - u(t)\\ p x dt = 0. (1) 
The following proposition is a consequence of Theorem [l] 

Proposition 2. Let 1 < p,q < 00. Let X\ C C X3 be real Banach spaces 
such that X\ is reflexive, the embedding X\ C X2 is compact and the embedding 
X 2 C X3 is continuous. If a subset Q C M p,q (I; X\, X 3 ) is bounded, then it is 
relatively compact in L P (I;X2). 
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rt2 






1 u(t) dt 






Jt! 





Proof. We apply Theorem [T] with X = X 2 . Let us hx < ti < t 2 < T and let 
v e G(ti,t 2 ). For u € Q we have 



(2) 



Thus G{t\,t*t) is bounded in X\ and therefore relatively compact in X 2 . 
It suffices to show the strong integral equicontinuity of Q. Let sup ug g II u IIlp(/-Xi) 
Af. We will use the Ehrling Lemma (see for instance [33] > Lemma 7.6). Let us 
fix e > 0. There exists C > such that for v e Ii we have \\v\\ x < graflMIx + 
C||w||^- . In particular, fixing ft, G (0. T),forn€5 and almost every t € (0, T— ft) 
we have |Ki + ft) - u(t)||» a < g^||u(t + ft) - u(t)\\ p Xi + C\\u(t + ft) - u (t)||* 
Integrating this inequality we get 



< 



< 



T-h 

e 

2PM 

e 



\\u(t + h) - u(t)\\ p x dt< 



T-h 



T-h 



\\u(t + h) - u(t)\\ p x dt + C \\u(t + h) - u(t)\\ p x dt< 



2M J 



T-h 



\\u(t + h)f Xi + \\u(t)\\ p x ^ dt + C I \\u(t + ft) - u(t)|&, dt < 



T-h 



\X 3 



T-h 



<e + C \\u(t + h) - u(t)\\ p x dt. 
Jo 3 

Now let 

sup u gcj Il u llsy9(7 jf3) — *~* ' ^£ P — Qi then by the Holder inequality, we 

have 



(3) 



T-h 



|K* + ft) - u(t)\\ p dt<T x -^ 



T-h 



IK* + ft) -«(*)!&, dt] . (4) 



If in turn q < p, then 



||u(i + /i)-u(t)|& dt < S« 



£-1 



T-h 



\\u(t + h) -u(t)\\ q x dt. (5) 



We estimate the last term in Q and ([5| from above (taking, if necessary, u(t) 
u(T), if t > T) 



f \\u(t + h)-u(t)\\ X3 dt< [ 
Jo i=Q Jil 

[T/h-2] h 

= V / \\u(t + ih+h)-u(t + ih)\\ 9 X3 dt = 

= Y] \\u(t + ih + h) - u(t + ih)\\ q X3 dt < Sh. 

Jo ,_ n 



[T//.-21 ih+h 



||ti(t + ^)-tt(*)||t dt 



(6) 
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Thus the last term in ([3]) tends to uniformly in u as h — > and, since e was 
arbitrary, we get the thesis. □ 



Remark 1. Note, that Theorem 3.2 in pQ is a consequence of above theorem. 
Compare also the Corollary 7.9 in (33] where the case p = 1 is excluded and X$ 
is assumed to have a predual space. 



4 Problem formulation and assumptions 

Let V C H C V* be an evolution triple, where V is a reflexive and separable 
Banach space and H is a separable Hilbert space with the embeddings being 
continuous, dense and compact. Embedding between V and H will be denoted 
by i. Furthermore let U be a reflexive Banach space on which the multivalued 
term will be defined. We use the notation V = L 2 (0,T;V), H = L 2 (0,T;H), 
U = L 2 (0, T; U) and W = {u £ V, v! £ V*}, where the derivative is understood 
in the sense of distibutions. Duality parings and norms for all the spaces will 
be denoted by the appropriate subscripts, for the space V no subscript will be 
used. Scalar product in H will be denoted by (•, •) and norm in K n by | • |. We 
consider the operator A : V — > V* and the functional J : U — > E such that the 
following assumptions hold 

H(A): (i) A is pseudomonotone, 

(m) A satisfies the growth condition < a + &||v|| for every v £ V 

with a > 0, b > 0, 

(in) A is coercive (A(v),v) > a||w|| 2 — /3|M|^ for every v £ V with a > 
and p > 0, 

H(J): (i) J is locally Lipschitz, 

(ii) dJ satisfies the growth condition \\£\\u* < C (1 + II M I|£/) f° r every u £ U 
and £ £ dJ(u) with c > 0. 

Moreover we assume that 

H : f £ V* and u Q £ H. 

We also impose the assumption concerning the space U 

H(U): There exists the linear, continuous and compact mapping i : V —¥ U such 
that the associated Nemytskii mapping I : M 2 ' 2 (0,T; V, V*) —>U defined 
by (lv)(t) = b(v(t)) is also compact. 

Finally we impose the last assumption 

H aU x' One of the following holds 

A) There exists a linear and continuous mapping p : H — > U such that 
for v £ V we have p(i(v)) = l(v). 

B) The constants a and c satisfy the inequality a > c||t||£/y.^\. 
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C) For every u £ U we have J°(m; -u) < d(l + \\u\\%) with d > and 
1 < a < 2. 

The problem under consideration is as follows 

find u £ W such that u(0) = uo and for a.e. t £ (0, T) we have 
«'(t) + Au(t) + L*dJ(iu{t)) 3 f(t). (7) 

The last inclusion is understood in the following sense 

there exists 77 £ V* such that u'(t) + Au(t) + 77(f) = /(f) for a.e. t £ (0,T) 
and (77(f), v) £ (dJ(Lu(t)), lv)u* x u for a.e. t £ (0, T) and v £ V. (8) 

Remark 2. The formulation ^ puts into a unified framework hemivariational 
inequalities originating from the initial and boundary value problem with mul- 
tivalued term defined on the problem domain (in this case we have multivalued 
source term, and U — H, see [221 HH [55] ) anc ^ 011 ^ e part °^ domain bound- 
ary d£l (this is the case if we have the multivalued, nonlinear and nonmonotonc 
boundary condition of Neumann-Robin type, U = L 2 (Tc) or U = i 2 (rc;K"), 
see [27] ). A detailed discussion as well as examples of problems which satisfy 
the assumptions will be given in Section [8] 

Remark 3. For the sake of simplicity of further argument the assumptions 
given above are not the most general ones under which the results hold. Possible 
generalizations include: 

• The dcpendance of A and J on time variable. Time dependent operator 
A for parabolic HVI is considered in [24] and the case of both A and J 
depending on time is considered in [27] (see Remark 8.21 in [33] on the 
Rothe method for the problem with the operator depending on time). 

• Instead of pseudomonotonicity one could assume that A is a sum of two 
operators, one of which is pseudomonotone and the second one is weakly 
continuous. Such weak continuity allows to take into account the nonlinear 
terms of lower order which are not of monotone type (see |13|). 

• More general coercivity conditions on A can be assumed. For instance 
(Av,v) > c\\v\\ 2 - a\\v\\ v - 7(t) with c > 0, a > and 7 £ ^(O^T) cf. 

HO- 

• The case when the space V is defined as L p (0, T; V) with 2 < p < 00 can be 
considered. Then we can assume more general growth conditions on A and 
J. For instance in [5B] it is assumed that ||A(t,t;)||v» < 0(t) + Ci\\v\\ p ~ l 
and that for r\ £ dj(x,£) we have \r}\ < c(l + IC^" 1 )- Note that J is 
defined typically as the integral functional J(u) = J j (x , u(x)) dx and 
assumptions on the integrand j are given. 

Remark 4. In this paper the abstract setting is considered. For a divergence 
differential operator of Leray - Lions type on a Sobolev space pseudomonotnicity 
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is implied by the appropriate Leray - Lions type conditions (see, for instance, |3J 
where conditions that guarantee pseudomonotonicity on W m,p {Vl), 1 < p < oo, 
m > 1 are considered). 

We conclude this section with the Lemma on pseudomonotonicity of Nemytskii 
operator with respect to the space M 2,2 (0, T; V, V*). Note that the proof of this 
lemma is analogous to the proof of Theorem 2 (b) in [3J (see also Proposition 
1 from [3U] and Lemma 8.8 in [33] for similar results). Lemma 8.8 of [33] is 
most similar to Lemma [l] but note that here no a priori bound in L°°(0, T; H) 
is needed and the assumption on the bound of 2-variation which is used here is 
weaker then the bound on 1-variation as in |33j . 

Lemma 1. Let A : V -> V* satisfy H(A) and let A : V V* be a Nemytskii 
operator for A defined by (Au)(t) — A(u(t)). Then if, for a uniformly bounded 
sequence {u n } C M 2,2 (0, T; V, V*) such that u n — > u weakly in V we have 
limsup„^ 00 (^lu„, u n — m)v«xv < 0, then Au n — > Au weakly in V* . 

Proof. It is enough to show that the thesis holds for a subsequence. By the 
generalized Lions Aubin Compactness Lemma (see Proposition [2| for a subse- 
quence (still denoted by n) we have u n — ¥ u strongly in %. Moreover, for yet 
another subsequence u n (t) — » u(t) strongly in H for a.e. t £ (0,T). We denote 
the set of measure zero on which the convergence does not hold by N. Now let 
us define £ n (t) = (Au n (t) , u n {t) — u{t)). We have 

> a\\u n {t)\\ 2 - (3\\u n (t)f H - \\u{t)\\(a + b\\u n (t)\\) > (9) 

> ?ll«nWII 2 - ^K(*)ll^ - a|K«)ll - ^\Ht)\\ 2 . 

2 let 

Now let C = {t 6 [0, T] : lim inf n _>oo (*) < 0}. This is the Lebesgue measur- 
able subset of [0, T]. Suppose that m(C) > 0, m being one dimensional Lebesgue 
measure. For every t £ C\N the sequence u n (t) has a subsequence (still denoted 
by n) which is bounded in V by ([9| such that lim n -n x {Au n (t), u n {t) — u(t)) < 0. 
Again for a subsequence we have u n (t) — > u(t) weakly in V, where the limit 
equals u(t) since we can consider only t (fc N. By the pseudomonotonicity of 
A we get < liminf„_ >00 (ylu„(i), u n (t) — u(t)}, which is a contradiction. So 
m(C) — 0, which means that liminf„_ i . 00 £ n (t) > a.e. on (0, T). From the 
Fatou Lemma we have 

P\\u\\ 2 H < [ liminf C„(t) dfc + y3||ul||r < / Iiminf(f„(t) + j8|K(*)||lr) <tt < 
Jo n ^°° Jo n ^°° 

< liminf / £ n (t) + f3\\u n {t)\\ 2 H dt < 



n— >-oo 







<liminf/ t n {t)db + P\\uf H <limsup / f n (t) dt + p\\uf H < p\\uf H . 

So Jq in{t)dt — y as 7i — y oo . Now note that |£„(t)| = Cn(*) + 2 C^(i) and 
^-(t) -> for a.e. t £ (0,T). Since, by for a.e. t £ (0,T) we have 
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Ut) + P\\u n (t)\\ 2 H > f(t) with / G Z^O.T), then^-(t)-p\\u n (t)f H < f~(t). 
Invoking Fatou Lemma again we have lim sup J £ - (t) dt < and further- 
more / T t;~(t)dt — > as n — > oo. We deduce that £„ — > in i 1 (0,T) and, 
for a subsequence (still denoted by the same subscript), £ n (t) — > for a.e. 
t £ (0,T). Since, for this subsequence, u n {t) — > u(t) weakly in V, then by 
pseudomonotonicity of A it follows that Au n (t) — > Au(t) weakly in V* and 
(Au n (t), u n {t)) — » (v4u(t), it(£)) . For any v £ V we have 

(.Ait, it - u)v*xv = / (Au(i),u(t) -v{t))dt = / lim (Ait„(i), u„(t) - u(i)) 

2 







/SIHIi + / lim ((Au n (t),u n (t) - v(t)) + P\\MWh) dt - 



n 



We can apply Fatou Lemma one last time to get 

(Au, it — v) y* x v < lim inf / (Au n (t) 7 u n (t) — v(t)) dt 



n 



= lim mi ((Au n ,u n - w)y*xv + (Au n ,u- v) v , xV ) < 

n— ¥ oo 

< liminf {Au ni u — w) v . x v- (10) 
Since w is arbitrary we obtain the thesis. □ 



5 The Rothe problem 

In this section we will work with a sequence of time-steps r„ — > such that each 
time step r„ > and the value T/r n is an integer, which we denote by N n . The 
subscipt n will be omitted in the sequel in order to simplify the notation, so we 
will write N, r instead of N n , r n . 

We define the piecewise constant approximation of the function / G V* . For this 
purpose we take the sequence of positive numbers e(r) — > and the sequence of 
mollificrs p e : K — > R which belong to C°°(E) and are nonnegative, supported 
on [— e, e] and f R p e (x)dx = 1. The function / is regularized according to the 
formula 

T / T -It 



u 



f e (t) = I Pe[t + e^^-s)f(s)ds. 



Note that f e £ C^O, T; F*) (see [33] , Lemma 7.2). The piecewise constant 
approximation for / is given by 

f T {t) := ! k T = fe(T)(kr) for te((i- l)r, fcr], fe G {1, . . . , N}. 

Following [33], Lemma 8.7, we have / T — > / in V* when r — > 0. Note (see 
Remark 8.15 in [33]) that the smoothing of / is not the only possible approach 
here. It is also possible to take the Clement zero-order quasi interpolant /™ = 
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We approximate the initial condition by elements of V. Let {mot} C V be 
a sequence such that uq t — > uq strongly in H and ||itoi-|| < f° r some 

constant C > 0. 

We define the following Rothe problem 

find the sequence {u^.}^ =0 C V such that u° = uor and 

( ^~ T ttT " 1 '«) +(^,«> + (3JK),i^'xt f 9 (/?,«) (11) 
for all w e y and fc = 1, . . . , JV. 

The above formula is known as the implicit or backward Euler scheme. Existence 
of solutions to the Rothe problem follows from the following 

Lemma 2. Under assumptions H(A),H(J),Hq, H(U) and H aux there exists 
tq > such that the problem Uj\) has a solution for t £ (0, to). 



Proof. We show that, given 1 G V, we can find £ V such that (111 
holds. We need to show that the range of multifunction V 3 v —> Lv = LJ:2L + 
Av + L*dJ(iv) e 2 V constitutes the whole space V* . We will use the surjectivity 
theorem for pseudomonotone operators (see for instance Theorem 1.3.70 in [T2]). 

We need to show that L is coercive (in the sense that lim||. u ||_ > . 00 mt "*^|^'j^" '"^ = 
oo ) and pseudomonotone. 

Claim 1. L is pseudomonotone. We verify this condition for all components of 
L separately. For this purpose we use Proposition 111 The operator ^ satisfies 
the conditions (£) — (Hi) trivially. As for i*dJ(m) the condition (i) follows 
from the fact that the Clarke subdiffercntial has nonempty, convex and (for 
reflexive space) weakly compact values. The condition (ii) follows from the 
growth assumption on dJ . In order to verify (Hi) let us take v n — > v weakly in 
V and £ n — > £ weakly in V* with € i*dJ(iv n ). Obviously iv n — ¥ lv strongly in 
U. Define r\ n G dJ(bv n ) such that £ n = i*T) n . By the growth condition H(J)(ii) 
it follows that, for a subsequence still denoted by the same subscript, rj n — > r\ 
weakly in U* . By the closedness of the graph of d J in U x topology (see 
[TO] . Proposition 2.1.5), we get r\ £ dJ(iv). Obviously £ = i*r\ and £ £ i*dJ(iv). 
Moreover = (iv n ,7] n )u* xU -> (w,ti)u*xU = (v,0> where by uniqueness 

convergence holds for the whole sequence. 

Claim 2. L is coercive. Assume that v* £ Lv. We estimate (v*,v) from below. 
For some rj £ dJ(iv) we have 

(v*,v) > -\\v\\h +a\\v\\ 2 -p\\v\\ 2 H + (r ] ,Lv) u * xU . (12) 

T 

We proceed for cases A),B),C) separately. For A) and B), by the growth 
condition 

(v*,v) >fl-p\ \\v\\ 2 H + a\\v\\ 2 - c(l + WlvWu^lvWu. 
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In the case A) we have \\iv\\1; < \\p\\% H m IMIh, so 

{V*,V) > (J^-P- C\\P\\C(H,U)] Nil + a \\ v \\ 2 - c lklU(V;I7)ll v l|- 

We require rn = . ,, ,} . In the case B) we get 

P + C I|P||C(H,[/) 7 ° 

(v*,v) > (J- - pj \\v\\% + (a- c||4||| (V;l0 )|M| a - c\\t\\ c(V;U) \\v\\. 

To have coercivity we need to set to = 4 . Finally if C) holds, then we get 

(v, i>v)u* xU > ~J°(lv; -lv) > 

> -d(i + uhi&) > -d - iitiiZ(v ;l oii«ir > -<*- f h 2 - c 

where C > depends on a,cr and Combining the last estimate with 

Q we get 

<«»> IHIIr + fNl a -d-C. 

Again setting r = 4 we get the desired property. □ 

Next lemma establishes the estimates which are satisfied by the solutions of 
Rothe problem. 

Lemma 3. Under assumptions H(A), H(J), H , H(U) and H aux there exists 
tq > such that for all r € (0, To) the solutions of Rothe problem (11) satisfy 



max \\vZ \\h < const 

fe=l,...,AT 

$^||u£-u£ _1 ||h < const, 
fe=l 
JV 

t^ ||u^|| 2 < const, 
fe=l 



(13) 
(14) 

(15) 



loii/i t/ie constants independent on t . 



Proof. We take v — u\ in ( 11 ), which gives for e > and k = 1, . . . , N 
i||«*||| + a||u*|| a + <^,t«*)tr.xir < 

T 

<^lklll, + ^ll/ T fc |l^+|lkl| 2 + ^(^- 1 ) ^) 

with ^ <= dJ{Lu k T ). We use the relation ||a|| 2 - (a, b) = ||a|| 2 /2 - ||b|| 2 /2 + \\a 
b|| 2 /2 to obtain 



2r 



P ) WKVh + - «rllfr + ( « - o ) IKIr + (16) 
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Recall that 

[-Cx-c\\ P \\l {Hm \\u k T f H -flKH 2 if A) holds, 

Lu k T )u, xU > I (-c\\i\\ 2 c{u v) - S)\\u k \\ 2 - C 2 if B) holds, 

l-C 3 - f ||^|| 2 if C) holds, 

where Ci > depends on c, a, ||£||jC(y ; [r)> 5 > is arbitrary, C2 > depends 
on c, S, \\l\\c(v;U) an d C3 > depends on d, a, cr, ||i>||,c(y ; m- From now on we 
proceed separately for the cases A), B) and C). In the case A) we take £ = f 
to get 



(17) 



- /3 - W)) II^Hh + - uj-ii, + 
+fll^ll 2 <^ll/ r fc ll^ + ^ll^- 1 ll 2 ff + c 1 . 

Summing above inequalities for fe = 1, . . . , n, where 1 < n < AT, we have 

n n 

IKIIl. + ^ll^-^ll^ + ^^ll^ll 2 ^ (18) 

fc=l fe=l 



< 2T Cl + 2r(/3 + cllpll 2 ^^) £ ll«*||2r + + Klllr- 

k=l 

Now if t < l/(4(/3 + c\\p\\j.^ HU - j j), by a discrete Gronwall inequality (see e.g. 

[33J (1.68)-(1.69)), we have Q-(JTIJ). 

In the case B), for S = e = (a — we get 

1 a\ IL.fe||2 i 1 ||„,fc „,fe— 1 1|2 , "~ C H i lli(t/,V) „, fe||2 



, r /3jii<n^ + -ii<-<- i ii^ + — ^f^n<r< (19) 

<C4||/ r fc ||^+^|h4- 1 ||l f + C5, 

where C4 = _ ,, h and C5 > depends on a, c, ||t|| %uy)- In analogy 



to the previous case we get ( 13 )-( 15 ] for t < 1//3. Bounds in the case C) are 
obtained in an analogous way. □ 



6 Convergence of the Rothe method 

We define piecewise linear and piecewise constant interpolants u T £ C([0, T]; V) 
and u T £ L°°(0, T; V) by the formulae 

u T (t) = f * - k+ l) u k T + (k- ^) u k - x for t E [(* - l)T,fcr], 
M T (i) = u k for a.e. i 6 ((fc — 1)t, fer]. 
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where k = 1, . . . , T/t. 

The sequences {u T7l }^ =1 and {u Tn }^ =1 are known as the Rothe sequences. 
Observe, that u T has a distributional derivative u' T G L°°(0,T; F) given by 

u' T (t) = Ut " t — for almost every t G ((fc — l)r, fcr). So, since solves the 
Rothe problem, we have for almost every t G (0, T) 

{u' T {t),v) H + (Au T {t),v) + (£, T {t),Lv) u , xU = {f T (t),v) for V£V, 

with u T (0) = u 0r and £ T (t) = ^ G 9J(«4) = dJ{iu r {t)) for < G ((fc - l)r, fer]. 
Defining the Nemytskii operator „4 : V — >• V* as (^4u)(t) = we have 



(u' T ,v)n + (Au T ,v) V * xV + (£t,w)w*xW = (/ T ,!))vxV for V G V. 



(20) 



Lemma 4. Under assumptions H(A), H( J), H , H(U) and H aux there exists 
To > suc/i that for all t G (0,to), f/ie piecewise constant and piecewise linear 
interpolants built on the solutions of the Rothe problem satisfy 



V < const, 

i«(0,T;H) < Const, 



c(o,T-,H) < const, 
V < const, 
v* < const, 
\Au T \\v* < const, 
|£r||w* < const, 
I«t||bv»(o,T;V) < const 

wzi/i i/ie constants independent on t. 

Proof. Estimates (|2~l])-(|2"3"|) follow directly from Lemma [3j since 



(21) 
(22) 
(23) 
(24) 
(25) 
(26) 
(27) 
(28) 



= T Ei=i 



,fe||2 



*r||i~(0,T;J?) — maXfc^.^jv || M T ll-tf 



and 



|m t ||c(o,T:H) < max fe=0: ...,7V \\K\\h- 



The simple calculation shows us that ||u T || v < T^tLn ll u rllv- 
with the fact, that ||u°|| < C/^/t, by Lem ma [3] gives (24). 



To prove ( 25 ) let us consider the inclusion ( 20 1 . We have 



This, together 



K||v* = SU P \{ u 't, v )u\ 

ll"llv<l 



sup 

Hlv<i 



(fr,v)v'xV ~ (Au T ,v) V . xV - I (S,r{t),iv(t)) U , xU dt 





< 



< WfAv + WJ \\Au T (t)\\ 2 v . dt+\\l\\ C(V;U) J UAtWu* dt < 



< UMlv. + j2a^T + 2bi\\u T \\ 2 v + \\i\\c(v,u) J^T + 2c^\\l\\1 {v \\u T ^) 
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Desired bound is obtained by (21). Estimates that appear in (29) prove also 
(26) and (27). It remains to prove (28). Let us assume that the seminoma 
BV 2 (0, T; V* ) of piecewise constant function u T is realized by some division 
= to < t\ < • ■ • < tk — T. Each tj is in some interval ((to., — 1)t, TOj-t], so 



u T (tj) — u T 3 with to = and toj. = N and m^+i > for i = 1, 
Thus ^ 



We use the inequality 



,JV-1. 



' 1 j _ „ m J-i 



< (rrij - rrij-x) ^ 



Thus 




(rrij - m,j-i) 



i=m,j — i-\-l 



' E i 

i=rrij — i-\-l 



< 



N 



N 



< ( £( TOj - m,_! - 1) ] £ K - «r _1 |lv' < ^E 



K(*)llv« 



The last term is bounded by ( 25 ) , which ends the proof. 



□ 

the problem 



Theorem 2. Under assumptions H(A), H(J),Hq, H(U) and H, 
£?J) has a solution u. Furthermore if u T and u T are piecewise constant and 
piecewise linear interpolants built on the solutions of the Rothe problem, then, 
for a subsequence, u T — > u weakly in W and weakly* in L°°(0, T; H) and u T — > u 
weakly in V and weakly* in L°°(0, T; H). 

Proof. From the bounds obtained in Lemma [4j possibly for a subsequence, we 
get 

u T — s- u weakly in V and weakly* in L°°(0, T; H), (30) 
u T — > u\ weakly in V and weakly* in L°° (0,T; H), (31) 
u' r — > W2 weakly in V*, (32) 
Au T — )■ 7? weakly in V*, (33) 
£ T -> £ weakly in W*. (34) 

A standard argument shows that =1*2- To show that u = u\ we observe that 



•>t\\v* 



kr 



N 

k-.l Ji llr 



(fcr - if 



eft = 



vllv* ' 



(35) 
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which means that u T — u T — > strongly in V* as r — > 0, and, in consequence 

U = U\. 

It follows that u T ^u strongly in £ 2 (0, T; H) and weakly in C([0, T];H). This 
also implies that Ug T = u T (0) — ¥ u(0) weakly in H, so u(0) = u a . 
A passage to the limit in |20| ) gives 

U + T) + L*£ = f. 

We observe that, by H(U), we have lu T — > lu strongly in U and, furthermore, 
for a subsequence iu T (t) — > iu{t) strongly in U for a.e. t <G (0,T). Moreover 
& — » £, weakly in L 1 (0,T;[/*). Since dJ : U — > 2 U has nonempty, closed 
and convex values and is upper semicontinuous from U furnished with strong 
topology into U* furnished with weak topology (see [IT], Proposition 5.6.10), 
by the Convergence Theorem of Aubin and Cellina (see [2J, Theorem 1, Section 
1.4), we deduce that € dJ(iu(t)) for a.e. t € (0,T). In order to show that 
u satisfies the inclusion ([7]), it suffices to prove that rj = An. To this end, let us 
estimate 

limsup(^u T , u T — u)v«xv < limsup(/ T ,u r — u)v*xv — (36) 

r^O r-s-0 

- liminf(?/ T ,M T - u) v * x v - liminf(^ T , l(u T - u))u»xu- 

t— >0 r— >0 



Since f T —t f strongly in V*, by (30), we get ]im T ^ ) . {f T ,u T — u)y*xv = 0. 
Moreover, since lu T — ¥ lu strongly in U, by (34), we have lim^o (£ T , l(u T — 
u))u*xu = 0. Now we observe that 

(u' T ,U T - 1t)v*xV = {u' T ,U T - M T )vxV + 

+ l -{\\u T (T) - u(T)\\ 2 H - K(0) - u(0)||^) + (u',u r - ujvxv, 

so, noting that (u' T ,u T — u T )v«xv > 0, we obtain 

lim inf (it' r , u r — w)y*xv > 0. 

Thus we have 

limsup(^4u r , u T — m)v»xv < 0. 
We are in a position to apply Lemma [TJ which gives r/ — Au. Thus u solves 

Remark 5. Note that we have also proved that any cluster point of u T and 
u T , in the sense (30)-(32), solves the problem ([T]). It is not known, however, 
whether there are solutions which are not limits of the interpolants built on the 
solutions of Rothe problem. 
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7 Uniqueness and strong convergence 

In this section we assume the strong monotonicity type relation for A and relaxed 
monotonicity on J. 

H(A)i : assumptions H(A) hold and A satisfies the monotonicity type relation 
(Au — Av,u — v) > mi||u — v\\ 2 — m2\\u — v\\ 2 H for every u,v S V with 
mi > and m 2 > 0, 

H{A)2 : assumptions H(A) hold and the Nemytskii mapping A : V — > V* is of 
class (5*+) with respect to the space M 2 ' 2 (0, T; V, V*), that is if u n — > it 
weakly in V and u n is bounded in M 2,2 (0, T; V, V*) then lim sup^^^ (Au n , u„— 
«)v«xv < implies that u n — > u strongly in V, 

H(J)i : assumptions H(J) hold and J satisfies the relaxed monotonicity con- 
dition (£ — 77, m — v)jj x u* > — m 3 \\u — for every «,d € F and £ € 
dJ(u),rj € dJ(v) with 7713 > 0, 

-ffconst : either if au2 . A) holds or mj > rn^\\L\\ 2 c( y .jjy 

Remark 6. The assumption H(A)i for the divergence differential Leray-Lions 
operator is guaranteed by appropriate Leray-Lions type conditions. For H{A)2 
to hold it suffices that the operator A is of class (S+), by an argument analogous 
to Theorem 2(c) in [3J. 

Remark 7. The relaxed monotonicity condition H(J)i (which is associated 
with the semiconvexity of the functional J) was already used to prove the 
uniqueness of solutions to the first order evolution parabolic hemivariational 
inequalities in [21] and second order ones in [28] . 

Remark 8. Note that H(A) 1 allows the case mi = 0, but if the inequality in 
H C onst holds, then it must be mi > 0. 

Theorem 3. Under assumptions H(A)\, H( J)i, Hq, H(U), H aux and H const , 
the solution to the problem is unique. 

Proof. Assume that Ui,« 2 are two distinct solutions to the problem Q. We 
have, for v g V and a.e. t g [0, T] 

(K - u 2 Y(t),v) + (A Ul (t) - Au 2 (t),v) + (£(t) - r]{t),Lv) UxU , = 0, (37) 

where £(t) € dJ(mi(t)) and n(t) g 9J(m 2 (t)) for a.e. t g (0,T). Taking 
v = U\ (t) — U2 (t) , we obtain 

~ "2(*)llff + (Aui(t) - Au 2 (t),ui(t) - u 2 (i)) + (38) 
+<£(*) - r]{t),Lu x {t) - Lu 2 (t))) UxU * = 0. 
Application of iJ(^4)i and H{J)x gives for a.e. £ g (0,T) 

^IM*) - u ^)\\h + ™i|M*) - «2MH 2 - (39) 
-m 2 |K(i) - ua(t)||2r - ms||t(«i(t) - u 2 (t))||^ < 0. 
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By Hconst we have either 
1 d 



. >f Ju 1 (t)-u 2 (t)\\ 2 H <m 2 \\u 1 (t)-u 2 (t)\\ 2 H , (40) 
or, in the case of H aux A), 

\j t \\ U ^) ~ U ^)W 2 H < (™2 + \\p\U(H;U))\\Mt) - Mt)\\ 2 H, (41) 

which, by the Gronwall lemma, gives the thesis. □ 

Remark 9. Under assumptions of Theorem [3j the convergences in Theorem [2] 
hold for the whole sequences u T and u T . 

Theorem 4. Let assumptions H(A)i, H(J), H^, H(U), H aux hold and the 
subsequences u T ,u T converge in the sense \3(fy -(32). Then u T — > u strongly in 
C([0,T]; H). If instead of H(A)\ we assume H{A) 2 , then u T — » u strongly in 
V. 

Proof. Let u T and u T be interpolants built on the solutions of the Rothe problem 
and let u be the solution to ^ obtained in Theorem [2j For v e V and a.e. 
t e (0, T) we get 

«(<) - u'(t),v) + (Au T {t) - Au(t),v) + (Ut) - V(t), lv}u, xU = (f T (t) -f(t),v), 

(42) 

where £ r (i) € dJ(iu T (t)) and r](t) € dJ(m(t)) for a.e. t £ (0, T). Choosing 
v = u T (t) — u(t), we get 

«(t) - u'(t), U T (i) - U T (t)) + ^IM*) - «(*)llff + 

+ (Azt T (t) - Au(t), u r (t) - u(t)) < 

< (f T (t) - f(t),u T (t) - u(t)> + ||m r (t) - M*(t)|M|&-(*) - 

Since (vf T (t),u T (t) - u T (t)) = \\^u' T (t)\\%(kT - t) > for any t € ((k - l)r, fcr) 
for a.e. i € (0,? 1 ) we have 

^IK(i) ~ «(*)lllr + (Au T (t) - Au(t),u T (t) - u(t)) < 
< \\iti T (t) - Lu(t)\\uU T [t) - rj{t)\\u, + 

+ </r(*) - f(t),U T (t) - «(*)> + («'(*), «r(t) - «r(*)>. (43) 

Using i?(^4)i and integrating the last inequality, for i € [O,? 1 ], we get 
1, 



2 : ||wt(*) - u(t)\\ H < m 2 J \\u T (s) -u(s)\\ H dt + 
+V / 2c||r(u r - u)\\u(2T + \\lu T \\u + \\iu\\u) + 

+ ll/r - /I|v*(||Wt||v + IMIv) + (u',U T - M T )v*xV + ^\\uq t - Uo\\%- 
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The Gronwall lemma gives the strong convergence of u T to u in C([0, T]; H). 
In order to obtain the strong convergence in V, let us integrate (43) over (0, T). 
We have 

-\\u T (T) - u(T)\\% + (Au T - Au,u T - !l)vxV < 
< V2c\\l(u t - u)|| M (2T+ \\lu T \\ u + + 

+ ll/r - /l|v*(ll S r||v + ||«||v) + (u',U T - U r )v>xV + dl^Or - "o||ff- 

Passing to the limit, we get 

limsup(.Ai2 T — Au, u T — u)vxv < 0. 

T-i-0 

The thesis is implied by H(A)2- □ 

Remark 10. If, in addition to assumptions of the Theorem |1J also H(J)\ 
and H const hold, then, by Theorem [3j the whole sequences u T and u T converge 
strongly in C([0,T}; H) and V respectively. 



8 Examples 

In this section we provide examples of that problem setup which are particular 
case of the general problem considered previously. Moreover, we present a simple 
numerical example. 

Problem settings We assume that f2 C K™ is an open and bounded domain 
with smooth boundary. The space V is either _ff 1 (51;]R' n ) with m £ N (possi- 
bly, but not necessarily, m = n) or its closed subspace (which originates from 
homogeneous Dirichlet boundary condition on C dfl). Furthermore let 
H = L 2 (tt;M. m ). Then the embedding i : V — > H is continuous and compact. 
We consider two examples. 

• Multivalued term is defined on f2. We specify A C to be an open 
subset on nonzero measure and fix d € N. Furthermore we assume 
that M £ L°° (A;£(]R m ;R d )). Now U = L 2 (A; R d ). The mapping t 
is defined by (lv)(x) = M(x)((w)\a(x)). We observe that l : V — > U 
is linear, continuous and compact. By Proposition [2] the embedding 
M 2 ' 2 {0,T;V,V*) C L 2 (0,T;H) is compact, which implies H(U). Defin- 
ing p : H — » U by (pw)(a:) = M(x)(v\\(x)), we see that A) of iJ aua; 
is satisfied. The solution exists under assumptions H(A),H(J) and H 
(Theorem^. Additional assumptions H(A)i and H(J)i imply uniqueness 
of solution by Theorem [3] and strong convergence of the Rothe sequence 
u T in C([0,T]; H). Furthermore, if H(A)2 holds, then the sequence u T 
converges strongly in V. 

As the special case we can consider A = fi, m = n = d and A4(x) = I 
(identity) for all x E il. Then we recover U — H, which gives existence 
results in spirit of |26j . 
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• Multivalued term is defined on the boundary of tt. We specify Tc C dfl 
disjoint with T D . We take Z = H s (Sl;M. m ) with S £ [|,1). The con- 
tinuous and compact embedding V — > Z is denoted by i and the trace 
operator is given by 7 : Z — >• L 2 (rc; R m ). Furthermore let d £ N and 
M £ L°°{T c ;C{R m ;R d )). Now [7 = L 2 {T c ;R d ). The mapping t is de- 
fined by (iv)(x) = M(x)(( i yiv)(x)). The mapping 1 : V — >■ U is linear, 
continuous and compact. The spaces V C Z C V* satisfy the assump- 
tions of Proposition |jj so M 2 ' 2 (0,T;V,V*) is embedded in L 2 (0,T;Z) 
compactly. Therefore the assumption H(U) is satisfied. Since claim A) of 
H aux does not hold in this case, in order to obtain the existence of solu- 
tions (Theorem[2|) we need to assume H(A), H(J), Hq and either B) or C) 
of H aux . Furthermore, if H{A)\ and H{J)\ hold, then a subsequence of 
the Rothe sequence u T converges strongly in C([0, T];H) (Theorem [4]) . If 
moreover H(A)2 holds, then we also have the strong convergence of u T in 
V. If furthermore the relation between mi and 7713 given by H const holds, 
then the solution is unique (Theorem [3]) and the whole Rothe sequences 
u T and u T converge strongly in C([0,T]; H) and V respectively. 

In the case m = d = 1 and J\4(x) = I we recover the results of [2"7] . 
If m = n > 1 and is the unit outer normal versor on the boundary 
dQ,, then two special cases are d = 1, A4(x)(a) — v{x) ■ a and d = m, 
M{x)(a) — a — (v(x) ■ a)v{x). We recover the cases of the boundary 
conditions given in normal and tangent directions, respectively. 

Numerical example. Let us take £1 = (0, 1). The problem under consideration 
will be 

t) for (x,t) £ n x (0,T), (44) 
u(0,t) = for t £ (0,T), (45) 
u x (l,t)£-dj(u(l,t)) for t£(0,T), (46) 
u(x,0) =u (x) for x£(0,l). (47) 

We set V = {v £ ^(0,1) : u(0) = 0} and i? = i 2 (0, 1). Taking t fc = kAt, 
u(x,tk) := u k (x) and u G V, the Rothe problem has the following form 

u W - u dx + / u k+i^ Vx ^ dx + dj(u k+1 (l))v(l) 3 0, 



At 

u°( x ) = Uo (a;). (48) 

The problem in each time step will be solved by the Galerkin scheme. Let V n be 
a subspace of V consisting of piecewise linear functions constructed on a uniform 
mesh xq — 0, . . . , Xi — iAx, ... ,x n = nAx = 1 in (0, 1) such that dim^ = n. 
Let furthermore u k (x) be approximated by X)"=i a i v ii where {vi}" =1 forms the 
base of V n given by the duality condition Vj(Xi) — 6{j. We assume that a k is 
the value of the solution at the last mesh point (x = 1). The values a k+1 satisfy, 
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for j = 1, ... ,n 



™ 1 f 1 
~J2 a i&t J Vi ^ v ^ dx + dj(a* +1 )v. 



Vi(x)vj(x) dx + / vKxtyJx) dx 



(1) 3 0. 



Calculating the integrals and denoting d — -ffis, for j = 2, , , . , n — 1, we obtain 

1 



and for j = n, we have 
1 



fe+1 ' - d I = q?_, -+aj-+a 



a.- 2 , „fc 1 



rl 6 TUj 3 T " J+1 6' ( ' 49 ' ) 



a 



fc+i 



fc+l 



1 



<')+!^y(«£ +1 )3«£-i 



i 



, i 



Finally for the left, Dirichlct, boundary point we have 

2d I n 



fe+i / 2 , ../• ^l x -d)=c$* ' 



i'2 



1 



:\ ' "J ' "* V 6 J 6 
We consider two examples of the locally Lipschitz functionals j: 



(50) 



(51) 









0,5 




\ r 





3JM 




0,5 




r 



2,5 -0,5 



Fi gure 1 ! Examples of multifunctions used as nonmonotone and multivalued boundary 
conditions Col 



for r < 



for r < 1 



_ J l-(r-2) 3 



= i t for r e (0, 1) j 2 (r) = 

\ for r > 1 [ | for r > 2. 

Their subdifferentials in the sense of Clarke are given by 



for r e (1,2) 



for r < or r > 1 
9ji(r) = { t for r € (0, 1) 
[0, 1] for r = 1 



for r < 1 or r > 2 
^. /2 (r) = <( 2-r for r e (1,2) 
[0, 1] for r = 1. 



(52) 



The graphs of <9ji and <9j2 are presented in Figure [T] Both potentials satisfy 
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Figure 2: Simulation for the potential j%. Plot of unique obtained solution is drawn. 



H(J). The potential j\ does not satisfy H(J)\ since its subdifferential has a 
nonmonotone jump. The potential ji satisfies H(J)\ since its subdifferential 
has a monotone jump and nonmonotonicity is Lipschitz. In the case of ji 
the question of multiplicity of solutions remains an open problem (however the 
numerical simulation below show that it is more likely that there are multiple 
solutions) and in the case of ji, a single solution is expected (at least as long as 
the inequality in H const holds). 

In both cases we take uq(x) = 2. The following scheme is used to find solutions 
of (49 1-(51 ). In every time step at most three solutions can be found: 



Assume that the element of dj(a^ +1 ) for which (50) holds is equal to 
i.e. we fall on the horizontal line in the graph of dj. Solve the system 
of n equations ( 49 )-{ 51 1 and verify whether obtained a^ +1 falls in the 
corresponding interval. 



Assume that the element of dj(a^ l +1 ) for which ( 50 1 holds is on the oblique 
line in the graph of dj. Solve the system of n equations ( 49 )-( 51 ) and verify 
whether obtained a^j +1 is in the corresponding interval. 

Assume that we fall on the vertical line in the graph of dj, i.e. ck„ 1 = 1. 
Solve the system of n— 1 equations (49) and (51). Then verify if dj(a^ +1 ) 
calculated from (50) falls in the interval [0, 1]. 



The simulations were run for At = 0.01 and Ax = 0.01. For the case of ji only 
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Figure 3: Simulation for the potential j\. Plots of two (respectively maximal and minimal 
one) of obtained many solutions are drawn. 

one numerical solution was obtained (i.e. in every time step only one of above 
three cases occurred). The result is presented in Figure [2] For the case of j\ 
many solutions were obtained (i.e. there were time steps in which more then one 
of above three cases occurred) . Figure [3] shows two solutions with respectively 
maximum and minimum value of chosen in each time step in which multiple 
solutions were found. 
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